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Wave functions of plane polygonal billiards are investi- 
gated. It is demonstrated that they have clear structures 
(superscars) related with families of classical periodic orbits 
which do not disappear at large energy. 

PACS numbers: 05.45.Mt, 03.65.Sq, 42.25.Fx 

A central problem of quantum chaos is an adequate 
description of different types of quantum systems which 
do not permit exact solution. For example, it is well 
accepted that eigenenergies of chaotic systems are dis- 
tributed as eigenvalues of random matrix ensembles [1] 
and their eigenfunctions are described by Gaussian ran- 
dom functions [2] . But much less is known when a model 
is neither chaotic or integrable. 

Particular intriguing examples are plane polygonal bil- 
liards (PPB) whose classical mechanics is surprisingly 
rich. When all their angles are rational with tt these mod- 
els are called pseudo-integrable (PI) because their clas- 
sical trajectories cover two-dimensional surfaces of genus 
g > 1 (see e.g. [3]). It was established numerically [4] that 
spectral statistics of PI models in many aspects resembles 
statistics of the Anderson model at the metal-insulator 
transition [5]. In particular, the nearest-neighbor distri- 
bution displays a linear level repulsion at small distances 
and an exponential decay at large separations. 

The purpose of this letter is to investigate wave func- 
tions of certain PI systems. It is found that they have 
superscaring property, namely, many of them have clear 
structures connected with families of classical periodic 
orbits which do not disappear at large energy. 

The scar phenomenon in chaotic systems is not new. 
The existence of structures near unstable periodic orbits 
in chaotic wave functions was established in [6] and later 
many works were done to clarify the subject (see e.g. 
[7] and references therein). PI systems differ in many 
aspects from chaotic and integrable systems and to the 
authors knowledge no conjecture about their wave func- 
tions exists in the literature. 

The main difficulty with analytical treatment of PI 
models is the strong diffraction on billiard corners with 
angles 7^ tt/u with integer n. When a bunch of par- 
allel classical trajectories hits these singular corners it 
splits discontinuously into two different bunches whose 
boundaries are called optical boundaries. Quantum me- 
chanics smoothes such singularities and associates with 
them scattering amplitudes [8] which due to discontinu- 
ous splitting of classical trajectories have different asymp- 



totics at large distances in different regions bounded by 
optical boundaries. Especially complicated case corre- 
sponds to multiple scattering when optical boundaries 
for different centers are close to each other. Up to now, 
due to singular character of such diffraction, it has only 
been proved analytically that the two-point correlation 
form factor for certain PI models takes at the origin a 
finite value different from standard statistics [9]. 

For PI systems classical periodic orbits form continu- 
ous families of parallel trajectories restricted by singular 
corners. Hence, waves traveling in such periodic orbit 
channels (POC) are influenced by infinite periodic ar- 
rays of singular diffr active centers. Fortunately, for the 
scattering on a staggered periodic array of half-planes 
this problem has an exact solution found in [10] by the 
Wiener-Hopf method. In [11] this solution was analyzed 
in the semiclassical limit of large energy and it was found 
that in the most singular case when the incidence angle 
(with respect to a plane formed by half-plane ends) is 
going to zero, all transmission and reflection coefficients 
tends also to zero except the 'elastic' reflection coeffi- 
cient (corresponding to the specular reflection from this 
plane) which goes to —1. It means that the fictitious 
scattering plane passing through singular diffraction cor- 
ners plays the role of perfect mirror with the Dirichlet 
boundary conditions. This mirror does not really exist 
but the multiple scattering on infinite number of parallel 
half-planes is equivalent to the reflection on this mirror 
plus corrections given by complicated formulas (see [11]) 
and governed by the perturbation parameter 



(1) 



where cp is the incidence angle with respect to the scat- 
tering plane, / is the distance between singular corners 
along the scattering plane, and k = \fE is wave momen- 
tum. Hence, when u ^ ^ (and /c ^ 00) the dominant 
approximation to the discussed multiple scattering prob- 
lem consists in treating scattering planes as true mirrors 
on which the total wave tends to zero. 

After unfolding each periodic orbit family in PPB cor- 
responds to an infinite POC restricted from both sides by 
straight lines passing through singular corners called sin- 
gular diagonals (SD). When a wave with a small u moves 
inside such channel, it reflects back and forth from SD 
as from perfect mirrors forming a propagating wave with 
zero boundary conditions on SD analogous to the Bor- 
rmann effect for scattering in crystals [12]. Therefore, in 
each POC one can construct the following quasi- states 
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called unfolded scar states which (i) obey the Dirichlet 
boundary conditions on SD and (ii) are periodic (or anti- 
periodic) along POC 

^^T^ sill (jm^ + sin {^riv) x{rt) ■ (2) 

Here ^ and r] {0 < r] < w) are coordinates, respectively, 
along and perpendicular to POC, / is the length of POC 
equal to the length of primitive periodic trajectory, w is 
the channel width, m, n > 1 are integers, and is a phase 
related with the choice of coordinates, x(^) (2) is the 
characteristic function of POC (x(^) = 1 or when x is, 
respectively, inside or outside POC) introduced to stress 
that scars states exist only inside POC. The bulk energy 
of such scar state is 

Em,n = [j) + ^-j . (3) 

Folding back the scar state (2) leads to a complicated 
expression ^f^^\x^y) which can be represented in a 
suitable expansion basis for any system. Folded scar 
states (i) obey the correct boundary conditions on bil- 
liard boundaries and (ii) fulfill the equation 
(A + Em^n)^ni,n\x^ v) = ^ everywhere except on SD. 

The above scar states exist if the perturbation param- 
eter (1) is small. As (/? ~ 7rn/wk the criterion of existence 
of a strong scar state with energy (3) is 

1 <n< rimax and rimax ^ wy/kjl . (4) 

This inequality implies that any POC in PPB will sup- 
port scar states with fixed n when ^ oo in marked con- 
trast with the scaring on unstable periodic orbits where 
contributions from individual orbits tend to zero in semi- 
classical limit. 

The total number of scar states depends on system 
considered. Only for special type of PI models called 
Veech polygons [13] (see also [1]) analytical calculations 
are possible. For such systems (i) the number of POC 
with / < L has the quadratic asymptotics 

N{1 < L) CL'/A (5) 

where A is the billiard area, C is a system dependant 
constant, and (ii) the width of POC with length / is 

w = -fA/l (6) 

where 7 < 1 is a constant taken from a finite set. 

Using (4) and these formulas one concludes that (i) 
scar states exist only for POC whose length is restricted 

I < /max with /max ^ k^^^ (7) 

(other channels are closed and cannot support propagat- 
ing waves with n > 1) and (ii) the averaged density of 



scar states ps is of the same order as the mean total den- 
sity of states p 

Ps = ^ S{E - Em,n) ^ ^rims^^^p. (8) 

scars ^<^max 

These results show that for PPB where (5) and (6) hold 
scar states are a good zero-order approximation to wave 
functions. As an illustration consider e.g. the right tri- 
angle with angle n/S. Its simplest POC corresponds to 
orbits perpendicular to the shortest side of the triangle. 
After unfolding it has a rectangular shape indicated in 
Fig. la. The folded scar state for this POC is shown 
schematically in Fig. lb. Dashed lines in this figure in- 
dicate maxima of the scar state. They have complicated 
form except at the right corner of the triangle where they 
form horizontal lines. 

In Figs. Ic-ld three true eigenfunctions of this triangu- 
lar billiard (with area 47r) are presented. The eigenfunc- 
tions were chosen because their energies are close to the 
scar energy (3) calculated with / = a and w = b where a 
and b = atan7r/8 are sides of the triangle. The charac- 
teristic horizontal lines corresponding to the scar picture 
(as in Fig. lb) are clearly seen in all these eigenfunctions. 
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FIG. 1. a) Unfolded scar state for the simplest POC of the 
right triangle with angle tt/S. b) Schematic folding of this 
state. Dashed lines indicate its maxima. Three solid lines 
show a region near SD where unfolded scar function tends to 
zero. c)-e) Eigenfunctions with energy E close to the scar 
energy Em,n- c) E = 407 A, ^50,1 = 407.6. d) E = 1015.97, 
E79,i = 1016.12. e) E = 1968.97, £;iio,i = 1969.15. 

These pictures are just a few examples (amongst many 
others) of clear scar eigenfunctions observed in trian- 
gular billiards. Complicated folding of POCs in such 
models makes it difficult to visualize scar states associ- 
ated with longer trajectories. This goal can more easily 
be achieved in another PI model, called barrier billiard 
(BE), which in the simplest case consists of a rectangu- 
lar billiard with the Dirichlet boundary conditions on all 
sides except the half of one side on which the Neumann 
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boundary condition is imposed [14]. In this model POCs 
are the same as for integrable rectangular billiard and 
are specified by two coprime integers M, TV. The POC 
length is / = ^/{2aMy + {2bNy where a, b are sides of 
the rectangle. The usual POCs for rectangular billiard 
are splitted (and restricted) by images of the singular 
point. For odd M each POC is divided into two POCs of 
width w = 2ab/L Both channels can support scar states 
but one requires odd m and the other - even m. When 
M is even, POC remains unramified and has the width 
w = Aab/L 

In Figs. 2 and 3 a few examples of high-excited scar 




FIG. 2. BB eigenfunctions strongly influenced by 1:0 and 
2:1 scar states, a) Folded scar states for the 1:0 (top) and 
2:1 (bottom) periodic trajectories. Thick line shows the part 
of BB with the Neumann boundary condition. Dashed lines 
indicate maxima of the scar state. Three solid lines show 
regions around SD where unfolded scar states tend to zero, 
b) Eigenfunction with E = 10209.55. The 1:1 scar energy 
^85 1 = 10209.65. c) Eigenfunction with E = 10017.57. The 
2:1 scar energy E453,i = 10017.67. 




FIG. 3. The same as in Fig. 2 but for 1:1 scar states. 

a) Folded scar states for 1:1 POCs with even and odd m. 

b) Eigenfunction with E = 10041.41. The 1:1 scar energy 
^347 1 = 10041.87. c) Eigenfunction with E = 10099.58. The 
1:1 scar energy ^348,i = 10099.82. 



states for BB with {b/a)'^ = >/5 + 1 and ab = 47r are 
presented. Black and white regions in these figures cor- 
respond to positive and negative values of eigenfunctions 
which are small in regions with irregular nodal patterns. 
This nodal domain representation is quite sensitive be- 
cause even a week noise changes drastically regular nodal 
pictures. Nevertheless, these (and many other) pictures 
show high-quality scar structures for BB. Fixing a POC 
and increasing the energy we always find reasonably good 
picture of the corresponding scar state close to the scar 
energy (3). 

More quantitative description of scar states is achi- 
eved by computing the overlap of folded scar states 
^^m,n\x^y) with exact eigenfunctions ^Ex{^->y) 



, {x, y)dxdy 



(9) 



In computations we fix n and choose m from the condi- 
tion of minimum of \Ex — Em,n\ (when m is kept fixed 
only one peak appears). In Fig. 4a we plot \Cm,n{E)\^ 
versus E with 2000 < Ex < 4000 for the 1:1 scar state. 
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FIG. 4. a) Overlap of exact eigenfunctions of BB with the 
1:1 scar state with (from bottom to top) n = 1, 2, 3, 4 (graphs 
for different n are shifted up for clarity by n — 1 units), b) 
Local density (10) for this overlap (graphs with different n 
are shifted up by 5(n — 1) units). Dashed lines indicate the 
best fit in the Breit-Wigner form (11). 

Spikes in such figures confirm that practically near all 
scar energies (3) there exists true eigenfunctions which 
have strong contribution from the given scar state. The 
analysis of these spikes reveals that their local density 



(10) 



averaged over different m can well be approximated by 
the Breit-Wigner distribution (see Fig. 4b) 



Pn{E) 



27r[(i?-e„(i?))2+r2(i?)/4] 



(11) 



similar to the one observed in random band matrices with 
preferential basis [15]. For the 1:1 scar state the best fit 
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gives ^ 3.5n^/V^ which agrees quahtatively with an 
estimate which can be obtained from [11] that for BB the 
total width of a given scar state in the leading order is 

For Veech billiards the density of scar states is of the 
same order as the total density of states (cf. (8)) and one 
is lead to to the conjecture that their eigenfunctions can 
be represented as a sum over scar states 

= E CmAm'^tT^^^y) ■ (12) 

scars 

For chaotic systems ergodicity according to Shnirelman's 
theorem [16] forces scar states to be rare. PI models are 
not ergodic and no a-priori objections exist which prevent 
this type of superscaring. 

Important characteristic of PPB wave functions is a set 
of participation ratios (see e.g. [17], [18] and references 
therein) 

R,{E) = {Y,\AmAE)?'')-' (13) 

where Am,n{E) are coefficients of the expansion of wa- 
ve functions in a suitable basis normalized such that 
Smnl^^^^P = 1. PI systems, as all systems with in- 
termediate statistics, should have some fractal proper- 
ties [17] and it is naturally to expect that Rq{E) ^^^^ 
]^Dq{q-i) ^]^]^ fixed Dq called generalized fractal dimen- 
sions (see e.g. [17], [18]). 

Under the simplest assumption that for scar states the 
local density of \Cm,n{E)\^^ is proportional to p^{E) with 
the Breit-Wigner form (11) of /)n(^), (12) and (11) with 
the above estimate for F^ give Dq = .5 confirming fractal 
character of BB eigenfunctions in momentum represen- 
tation. In Fig. 5 we plot Rq for g = 2 and g = 3 for BB 
computed directly from the expansion of eigenfunctions 
into trigonometric series. The best fits R2 ^ 2.52V^ and 
Rs ^ A.Tk very well describe the data in the given inter- 
val of energy which means that for BB D2 ~ ~ .5 
in accordance with the above estimates. The assumed 
value of spectral compressibility for BB, % = .5, [14] is 
close to the spectral compressibility x(^2) numerically 
computed at the point D2 = for the critical power-law 
random band matrix model (cf. Fig. 2 of [18]). 

In summary, we have argued that strong diffraction in 
PPB leads to the formation of new type of long-lived res- 
onant states (scar states) propagating inside POC and 
refiecting from SD as from perfect mirrors. Many true 
eigenfunctions of PPB have surprisingly clear structures 
associated with such states even at high energies. It fol- 
lows from our results that PPB are the best models of 
scar phenomenon. For good PI models (Veech billiards) 
the density of scar states is a constant and they can be 
considered as the basis of perturbation expansion. A 
week residual interaction between them (neglected in this 
paper) forms true eigenstates and leads to intermediate 



character of spectral statistics for these models. It ap- 
pears that this interaction shares many features with crit- 
ical random band matrix model. We have also checked 
that BB wave functions have fractal properties in mo- 
mentum space. 
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FIG. 5. Participation ratios R2 (top) and R3 (bottom) ver- 
sus energy for BB. White lines indicate the fits R2 ~ 2.52v^ 
and R3 ^ 4:.7k. 
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